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Abstract We derive the expression for the ponderomotive force in the real 



C/3 

C^ ' multicomponent magnetospheric plasma containing heavy ions. The ponderomotive force 



considered includes the induced magnetic moment of all the species and arises due to 
inhomogeneity of the traveling low-frequency electromagnetic wave amplitude in the 



C/3 

^ . nonuniform medium. The nonlinear stationary force balance equation is obtained taking 



into account the gravitational and centrifugal forces for the plasma consisting of the 



>■ ! electrons, protons and heavy ions (He"'"). The background geomagnetic field is taken for the 

cn I dayside of the magnetosphere, where the magnetic field have magnetic "holes" (Antonova 



and Shabansky 1968). The balance equation is solved numerically to obtain the nonlinear 



O ■ density distribution of ions (H"*") in the presence of heavy ions (He"*"). It is shown that 

for frequencies less than the helium gyrofrequency at the equator the nonlinear plasma 

'k>( ■ density perturbations are peaked in the vicinity of the equator due to the action of the 

c^ ■ ponderomotive force. A comparison of the cases of the dipole and dayside magnetosphere 

is provided. It is obtained that the presence of heavy ions leads to decrease of the proton 
density modification. 

Keywords Earth ■ Magnetic fields ■ Plasmas ■ Ponderomotive force ■ Waves 
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1. Introduction 

A considerable attention lias been paid to study an influence of heavy ions (mainly 
helium and oxygen) on the generation and dynamics of electromagnetic ion cyclotron 
(EMIC) waves in the frequency range 0.5 to 5.0 Hz traveling along the geomagnetic fleld 
lines. Using the searchcoil magnetometer and particle spectrometer data on board of the 
GEOS 1 and ATS 6 satellites, Young et al. (1981), Mauk et al. (1981) and Fraser et al. 
(1982) came to the conclusion that EMIC waves are strongly controlled by the dynamics of 
heavy ions. Later on, Kozyra et al. (1984) and Fraser et al. (1992) obtained similar results, 
using the ISEE 1 and 2 data. 

The propagation of the EMIC waves along magnetic fleld lines in a heavy ion rich 
plasma is characterized by the reverse of polarization and the splitting of the wave spectrum 
into two branches, the high-frequency branch uj > Uch {^ is the wave frequency, Uch is the 
heavy ion gyrofrequency) and the low-frequency branch u < Uch- The two branches are 
separated by a stop-band. Observations on board of the GEOS 1 and 2 satellites (Young 
et al. 1981) have shown that the EMIC wave spectra are concentrated in the vicinity of 
the equatorial He"^ gyrofrequency. These observations also showed that there is an inverse 
connection between the increase of the He^ concentration and the appearance of the Pel 
events on the ground that conflrms the important role of the He^ ions for the generation 
and propagation of the EMIC waves. 

As is well-known, ponderomotive forces induced by the EMIC waves contribute to the 
plasma balance in the Earth's magnetosphere (e.g. Allan 1992; Guglielmi et al. 1993, 1995; 
Guglielmi and Pokhotelov 1994; Witt et al. 1995; Pokhotelov et al. 1996; Allan and Manuel 
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1996; Feygin et al. 1998; Nekrasov and Feygin 2011, 2012). In studies mentioned above, the 
magnetosphere has been assumed to contain only one ion species (H+). In this paper, we 
explore the effect of the ponderomotive force in the real multicomponent magnetospheric 
plasma. In the numerical analysis, the plasma consisting of the electrons, protons and heavy 
ions (He+) is considered. The background geomagnetic field is taken for the dayside of the 
magnetosphere, where the magnetic field have magnetic "holes" (Antonova and Shabansky 
1968). The stationary nonlinear balance equation is solved numerically to obtain the 
nonlinear density distribution of ions (H+) in the presence of heavy ions (He"*"). Without 
heavy ions, this problem has been treated by Nekrasov and Feygin (2012). 

The paper is organized as follows. In Sect. 2, we describe the geomagnetic field model 
of the dayside magnetosphere. The ponderomotive force in the multicomponent plasma is 
derived in Sect. 3. The stationary force balance equation is considered in Sect. 4. In Sect. 
5, the results of numerical calculations of the balance equation for the nonlinear plasma 
modification in the curvature geomagnetic field are represented. Conclusive remarks and 
discussion are given in Sect. 6. 



2. Geomagnetic field model of the dayside magnetosphere 

We here apply a model of the Earth's magnetic field given by Antonova and Shabansky 
(1968). This model has been shown to be in a good agreement with magnetic field 
observations by the HEOS 1 and 2 satellites in the dayside magnetosphere (Antonova et al. 
1983). According to this model, the geomagnetic field is described by a superposition of 
two dipoles: the internal dipole with the magnetic moment M and the additional external 
dipole with the magnetic moment kM {k is a constant parameter) disposed at the distance 
a (measured in the units of the Earth's radius Re) on the dayside of the magnetosphere 
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along the Earth-Sun hne from the position of the original dipole. The latter dipole models 
the distortion of the geomagnetic field caused by the solar wind pressure. The model by 
Antonova and Shabansky (1968) is sufficiently simple and convenient for an analysis of 
geophysical phenomena in the dayside magnetosphere. We have described this model in our 
paper Nekrasov and Feygin (2012). However, for convenience of reading, we also give the 
main moments here. 

In the spherical coordinate system, magnetic field components for the two-dipole model 
by Antonova and Shabansky (1968) in the meridional noon- midnight plane have the form 

Br = ^tt, (1) 



ipO 



B 



B^ = ^v^r^^/J, (2) 



r 



where x = sirup, ip is the geomagnetic latitude, r is measured in the units Re, and Be 

= 0.311 G is the equatorial magnetic field at the Earth's surface. Coefficients a and /3 are 

the following: 

a = 1 ^ s,^ , (3) 

2 (a2 + r2 - 2arv^r^^) 

kr^ [y/l - x^ (a^ + r^) - ar (2 + x^)] 

Vl — x'^ (^2 + r^ — 2arVl — x'^^ 

The total magnetic field B in an arbitrary point of the field line can be defined from (1) 

and (2) as 

B = ^ [AxW + (1 - x') /32] '/' . (5) 



An equation for the field line is determined by dr/rdip = Br/ B^ or 

dr 2xr a 

dx 1 — x"^ (3 

The magnetic field B described by (5) has along the near boundary field lines two 
minima located symmetrically relative to the equator (Antonova and Shabansky 1968). 



(6) 



-5- 

When a tends to infinity, the values a and /3 tend to 1 (see 3 and 4). In this case, we have 
a transition to the one-dipole approximation. The dependence of the geomagnetic field (5) 
on X for different distances from the Earth is shown in Fig. 1, where L is the Mcllwain 
parameter. 



3. Ponderomotive force in the multicomponent plasma 

3.1. Ponderomotive force due to the longitudinal inhomogeneity of the wave 

amplitude 

To find the ponderomotive force of electromagnetic waves in the multicomponent 
plasma, we use the results of the paper by Nekrasov and Feygin (2006). The equation for 
the nonlinear slow velocity (vj2) is given by 



dt rrii 



((E2) + (F,2)), (7) 



"J 

where Qj and rrij are the charge and mass of the species j, (E2) is the slow nonlinear electric 
field, (Fj2) is equal to 

in ■ 1 

(F,2) = --^ (v,i ■ Vv,i) + - (v,i X Bi) , (8) 

qj c 

the angle brackets () denote the time- averaging over fast oscillations. The subscript 1 in (8) 
relates to the linear perturbations of the velocity, v^i, and magnetic field, Bi. 

We are interested in the nonlinear motion of a plasma along the background magnetic 
field. We assume that the latter has the z-direction. In the case of the EMIC waves also 
traveling along the magnetic field, the equation (7) for the longitudinal velocity (f j2z) takes 
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dt rrijU^^^ 

and the value {Fj2z) defined by (8) is given by 

zrrij auj [auj — uJcj) oz 

where a; > is the wave frequency, a = ±1 denotes the left (+) or right (— ) polarization 
of the wave, Upj = [Anrijq^ /rrij) is the plasma frequency, rij is the background number 
density, Ucj = qjB/mjC is the cyclotron frequency, c is the speed of light in vacuum and 
El = Re Eio exp (i f^ kdz — itot^ is the wave electric field with the nonuniform amplitude. 
When obtaining (9) and (10), we have used (11)-(13), (16) and (17) in Nekrasov and 
Feygin (2006), have taken into account that Vjiz = for waves under consideration and 
have assumed that d/dz ^ kd/udt, where d/dt is applied to the wave amplitude. The 
right-hand side of (9) (without rrij) is the particular (see below) ponderomotive force acting 
on species j. 

Using (9), we can define the total particular ponderomotive force in the multicomponent 
plasma as 

Taking into account (10), the expression (11) gets the following: 

87r -^ oOcj [cruj — cOcj) oz 

where we have neglected small terms proportional to rrie, which describe the ponderomotive 
force acting on the electrons. 
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3.2. Ponderomotive force due to the transverse inhomogeneity of the wave 

amplitude. Magnetic moment 

It is known that if tlie medium lias the own magnetic moment (magnetic field) M and 
is embedded in the external nonuniform magnetic field B, then this medium is subjected by 
the action of the force 

Fm = (M-V)B (13) 

(see e.g. Landau and Lifshits 1982). The EMIC waves propagating along the geomagnetic 
field generate the nonlinear magnetic field that can be found from (3), (16) and (17) in 
(Nekrasov and Feygin 2006). Neglecting in (16) (Nekrasov and Feygin 2006) small terms 
and passing to the space-time representation, we find the equation for the transverse 
nonlinear electric field {E2x) 

vv,-^)(E,,)^^^(a>, (14) 

where V^ = d'^ jdy^ + d"^ jdz^, c^ /c\ = Ylj ^Ij/^lj^ ca is the Alfven velocity and 

d{Cx) _ c ^ <• 9(E2) 



Ett^^^^ (15) 



dt 2B^{auj-u,,y dy 

(see (17) in (Nekrasov and Feygin 2006)). 

We further assume that d'^/dy'^ 3> d'^/dz'^ and c\d'^ /dy"^ ^ d"^ /dt^. Then from (14) 



and (15), we obtain 



Ly <■ ^(E^) (16) 



d{E2.) _ 1 ^ <■ 9(E2) 



dy 2cB ^ (aco - u,jf dt 

From Faraday's equation 

Vx (E.) = -i^, (17) 

c ot 

we can find the induced nonlinear magnetic field (B2). Substituting (17) into (16), we have 
in our case for the multicomponent plasma 



(^-) = -^E7-^(E?>- 



2^ ~^ (o"^ - (^cjY 



2 

1^ 



J 



We note that this expression for the two-component plasma (electrons and ions) has been 
derived by Nekrasov and Feygin (2005). 

The magnetic moment M is connected with the induced magnetic field as 47rM = (B2). 
Using this relation, we obtain in our case the expression (13) for Fm along the magnetic 
field 

This expression for the two-component plasma has also been obtained by Nekrasov and 
Feygin (2005), but in a different way. However, we see here under which conditions (18) 
takes place. 



3.3. Total ponderomotive force in the nonuniform medium 



The total ponderomotive force Fpz is equal to sum of (12) and (19) 



pz 



-E 



Tp \ PJ 



Ucj (au - Ucj) 



d (E?) u^j dB 



dz [auj — Ucj) Bdz 



n) 



(20) 



where we have neglected the contribution of the electrons into (18). We further assume 
that the wave amplitude depends on the z-coordinate because of the longitudinal medium 
inhomogeneity. In this case, the amplitude of the wave Eiq in the WKB-approximation is 

— 1/2 

proportional to Nz , where Nz is the refractive index equal to 

Nl = l-y^^i^ r (21) 

^ Ucj [cruj - u,,j) 

(see e.g. (14) in (Nekrasov and Feygin 2006)). When obtaining (21), we have used the 
condition of quasineutrality ^ • qjTij = 0. We see that the first term in the square brackets 
in (20) can be written in the form 
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Calculating the value d(^\) /dz through dN'^/dz, using (21) and substituting the resuh 
into (20), we find 



{^D sr^ ^Ij (dlnrij 



^^ IGtt ^ Ucj {uJcj — crcd) \ dz 



'' ^^ — '" au 2uci 1 1 dlnB 

+ ^ 



UJcj — (JCJ (cucj — o"w) N^ 



dz 



(22) 



Thus, the ponderomotive force in the multi-component plasma is the sum of the 
ponderomotive forces for each ion species. 



4. Stationary force balance equation 

From equations of motion for the species j in the second approximation on the wave 
amplitude averaged over fast oscillations, we can obtain the force balance equation in the 
stationary state. The total equation of motion contains except the electromagnetic force 
also the gradient of thermal pressure, the gravity and centrifugal force (e.g. Lemaire 1989; 
Persoon et al. 2009). 

We assume for simplicity that all the species have the equal temperatures 
Tj = T = const. Consider the three-component plasma consisting of the electrons (e), 
protons (i) and heavy ions (h). In the stationary case, adding the corresponding equations 
of motion for each species, we obtain for the nonlinear time-averaging density perturbations 
nj2 the following force balance equation along the magnetic field line: 



TVii 



2ni2 + {1 + — ]nh2 



(mjni2 + mhnh2) {g\\+nn\ft^ Ret cos if) + Fpy, (23) 



where the subscript || denotes the local z-direction. Here, we have used the condition of 
quasineutrality and the equality qi = —qe- In (23), the value gw = g ■ b = —gEBr/r'^B is 
the longitudinal gravitational acceleration {gE = 9.8 m sec~^), Q is the Earth's rotation 
frequency and n^n = nf7 ■ b, where n^ and b are the unit vectors along the centrifugal force 
and magnetic field, respectively. The value n^n is equal to nj^y = (1 + (3 /2a) {Bj./B) comp 
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(see Nekrasov and Feygin 2012). We note that in the last paper in (14) cos^ y? stands 
erroneously instead of cos ^. 

In (22), we will express the value (E^) through the amplitude of the wave magnetic 
field at the equator i?io. From Faraday's equation, it is followed that Nl (E^) = (B^), 
where (E^) oc A''^-'^. Thus, we obtain that (E^) = B\q {Nzq^z)~ for circularly-polarized 
waves. Here and below, the subscript relates to the values at the equator. The operator 
V|| in (23) is defined by the relation Vy = b ■ V and has been found in (Nekrasov and 
Feygin 2012) 

V|| = ^Re'v-'^'-^, (24) 

where 

We further connect the perturbation nh2 with nj2- From the nonlinear continuity 
equation and equation of motion (9), we can obtain an estimation of this connection 

m, H {1- i^,of ) 

where P = pi/pio, H = ph/pio, pi = niirii, ph = rjihUh, i^m = au/udo, vm = cru/ucho- 

For the three-component plasma, substituting (22), (24) and (26) into (23), we obtain 
the following equation: 

? = T f ^1 - ^ V + T^2 (^3 + ^4 + ^5 + A) , (27) 

ax A \ ax J A 

where 6 = Pi2/pio, Pi2 = mini2 and 

V qj mh P {1 - u^of) 
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In (27), we have introduced the notations 



^ 2RExafigeff , . 

/i2 - 



As 

Aa 



1 dP 



Pl/2 dx ' 

(1-^^iof) 1 dH 
(1 - ^.of ) i^^/' dx ' 



where c^j = 2T/mj and 

^e// = ^s + 9cf = 9E-il + /3/2a) fi^P^r^ cos^ if, (30) 

We consider that A^f ^ 1. The value /iq is obtained from (31) at a; = 0. When obtaining 
Ai, we have used (1), (5), (6) and (25). 

For the equihbrium mass density of H+ (i) and He+ (h), we take the power law form 
to describe the longitudinal field line distribution, pj oc r~^, where j = i, h. For the large 
distances from the Earth's surface which we consider below, the choice 7 = 1 is the best 
one to be appropriate to experimental data (Denton et al. 2006). Thus, we set 

Pj (^) = pjo (1 -x^y . 

This formula can be applied up to (/? ~ ±50 — GO'' (Denton et al. 2006). In this case, 
P = (1 — x"^) and H = Hq {1 — x"^) . These values will be substituted to (28) and (29). 
We note that in our case ruh = -imi. 



-12- 



5. Numerical analysis 

Equation (27) for the two-component magnetospheric plasma has been analyzed in 
(Nekrasov and Feygin 2012) by the Runge-Kutta method, making use of the boundary 
condition 6 = 0.06 at a; = 0. Thus, we have obtained effects of the action of the 
ponderomotive, gravitational and centrifugal forces on the plasma density distribution at 
X ^ 0, at the same time perturbations at the equator remained unchanged. In this paper, 
we set the boundary condition 6 = at x = 0.8. This case permits us to consider the 
plasma density modification in the region of the equator. For numerical calculations, we 
have used the same parameters for the magnetic field model as that in (Nekrasov and 
Feygin 2012): a = 33 and A; = 13 (except for Figs. 3a,b). These parameters correspond 
to the dayside boundary of the magnetosphere at the distance IORe obtained by HEOS 
1 and 2 satellites (Antonova at al. 1983). Other parameters are the following: c^- = 10^ 
m^ sec~^, piQ = 1.67 x 10^^° kg m~^ for all L near the midday boundary of the Earth's 
magnetosphere, where the plasma density depends weakly on L (Chappel 1974; Carpenter 
and Anderson 1992). In all the numerical calculations, we have taken a = +1 and assumed 
that Uho < 1 to avoid the singularity 1 — Uh = 0. As heavy ions, we use He"*". 

Figures 2 depict the distribution of the normalized nonlinear proton density S = Pi2/pio 
along the field line for different values of Hq = 0, 0.5, 1.0 at L = 6, 10 for Bio = 10^^ G. We 
have set z/jo = 0.1 to have Uho < 1. We see that the plasma density perturbation increases 
in the direction of the equator due to the action of the ponderomotive force (without the 
latter see Fig. 5). The equatorial value of 6 decreases with increasing of Hq because of 
increasing of A and p with Hq (see (27)-(29) and (31)). 

Figures 3 represent a difference in the (5-distribution for the dipole (k = 0) and 
two-dipole {k = 13) geomagnetic field at L = 8 and Hq = 0, 1.0 for Biq = 10^^ G. We 
see that the peak value of 6 at the equator is smaller for A; = 13 than for k = 0. This 
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can be explained by increasing of the equatorial geomagnetic field B due to k ^ (see 
(3)-(5)), which stands in the denominator of A2 (see (29)) and decreases the ponderomotive 
force. The presence of heavy ions also results in a decreased influence of the ponderomotive 
force on the ions due to increase of the parameters A and fi (as in Fig. 2). Therefore, for 
Ho = 1.0, the values of 6 are smaller than for Hq = 0. We note that the curves for Hq = 0.5 
are here and in Figs. 4 situated between the curves with Hq = 0, 1.0 (except Fig.5). 

Figures 4 show the dependence of the relative density perturbation 6 on x for different 
values of the wave amplitude Biq = (1, 3, 5) x 10^^ G for L = 10 and Hq = 0, 1.0. It can 
be seen that increased wave-amplitude yields increased values of 6 along field lines. The 
smaller values of 6 at Hq = 1.0 are explained in the description of Figs. 2a, 2b and 3a, 3b. 

Figure 5 describes formally the case -Bio = to show the role of the gravitational and 
centrifugal forces (see (30)). It is followed from (27) that gravitational force increases the 
plasma density (not only perturbations!) with larger distances from the equator. This 
is shown by the curve 1. The centrifugal force has the opposite direction (see (30)) and 
drives a plasma to the equator. Far from equator, the gravitational force is larger than 
the centrifugal one (see also Sect. 6). In the vicinity of the equator, on the contrary, the 
centrifugal force is larger. This results in a peak of plasma density at the equator. The 
curve 2 shows the corresponding distribution. The equatorial values of 6 for the curves 1 
and 2 are decreasing with increasing Hq. The boundary condition for the curves 1 and 2 is 
6 = 0.01 at a; = 0.8 since to solve (27) for A2 = it is necessary to have a nonzero boundary 
condition for 6. 

All the figures demonstrate that the density perturbation is peaked at the equator due 
to ponderomotive and centrifugal forces. Such peaking is observed in the magnetosphere of 
the Earth (Denton et al. 2006). 
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6. Conclusion and Discussion 

We conclude by summarizing main results obtained in this paper: 

1). We have derived the general expression for the ponderomotive force induced 
by electromagnetic ion cyclotron waves in a multicomponent plasma containing different 
species of ions. 

2). The correct equation of the force balance along the magnetic field lines, which 
contains the perturbed thermal plasma pressure together with the ponderomotive force, 
having the same (second) order of magnitude, has been considered. 

3). We have investigated the effect of the ponderomotive force on the perturbed ion 
(proton) density distribution in the presence of helium ions. 

4). It has been shown that for frequencies less than the helium gyrofrequency at the 
equator the nonlinear plasma density perturbations are peaked in the vicinity of the equator 
due to the action of the ponderomotive force. The maximum of the ion (proton) density 
perturbation has been obtained to decrease with increasing of the heavy ion (He^) mass 
density. 

5). We have obtained that larger wave-amplitudes inducing more ponderomotive force 
result in larger 5-perturbations. 

Theoretical and numerical results given above in this section are the main points 
of our exploration. The derivation of expression for the ponderomotive force in the 
multicomponent plasma containing different species of heavy ions is particularly important. 
This permits us to consider the influence of wave perturbations on the redistrbution of 
magnetospheric plasma in the direction of the equator in conditions close to real ones. 
As we see, the equatorial value of the proton density distribution should drop for more 
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abundance of He"*". 

We have considered the relative role of the gravitational and centrifugal forces when 
the wave activity is absent, Biq = (Fig. 5). We see from (30) that the gravitational and 
centrifugal forces along the magnetic field lines have opposite directions. At some latitude 
iPeq, both forces are equal to each other. Taking into account (6), the estimation of ipeq for 
a ~ /3 ~ 1 is the following: 

where toRe is the equatorial distance of the given magnetic field line from the Earth's 
center (rg = L). At (/? > (peq-, the gravitational force is larger than the centrifugal one. 
When if < (/9eg, the centrifugal force is dominant. For tq = 8, we obtain cosy^eg ~ 0.89 or 
ifeq ~ 27°. When tq = 10, we have cos ipeq ~ 0.82 and ipeq ~ 35°. 

From (27), we can roughly estimate the relative contribution of the ponderomotive 
force Fp and of the gravitational and centrifugal forces -F^+c as 

^^ - ^^^^ (33) 



Fg + c Ai5 

Assuming that (5-perturbation is due to the wave action, we estimate in our case 5 ~ A2. 

Substitution of 6 in (33) gives 

Fp ^3 _ 4^0 



rNJ 



Fg + c ^1 RsdE 

For tq = 10 and 4 = 10^ m^ sec~^, we find Fp/ Fg+c ~ 160. We see that the ponderomotive 
force is dominant because Ai <^ 1 (for x < 1). We note that a change of the boundary 
condition by the nonzero valus of 5 in Fig. 2b does not influence on the form of the curves. 
Thus, we can expect substantial increases in (5-perturbations during more geomagnetically 
active periods (Figs. 4a, 4b). This result is qualitatively consistent with observational data, 
such as Denton et al. (2006), which show a larger equatorial mass density peak with larger 
wave amplitudes (see their Figure 12). As a possible reason for this effect, Denton et al. 
(2006) indicate a role of the ponderomotive force in driving ions up the field line. 
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We have investigated the nonhnear plasma density redistribution for the Antonova 
and Shabansky (1968) model of the geomagnetic field B given in Fig. 1. This choice is 
justified by its sufficient simplicity. In addition, the model by Antonova and Shabansky 
(1968) provides a complete analytical description of the geomagnetic field on the dayside 
of the Earth's magnetosphere and is in accordance with satellite measurements. The use 
of a more complex geometry for B such as T96 (Tsyganenko 1995) could be justified for 
taking into account some details (the local time, position, the solar wind parameters etc.). 
However, the results obtained will qualitatively remain the same because the main form of 
the magnetic field lines does not essentially change. 

In this paper, we have applied the theory of perturbations to obtain the ponderomotive 
force expression and the force balance equation. We assumed that nj2 ^ njQ. The value 
6 for parameters used is in the region 5 ~ 2 x 10^'^ ^ 6 x 10^^ (see Figs. 2-4). Thus, 
perturbations are small in comparison with the background density. During active wave 
periods, density perturbations can be of the same order of magnitude as background 
values (see Denton et al. 2006). Therefore, there is some problem to compare theoretical 
and experimental results. However, for the limiting estimation, we could take such wave 
amplitudes, for which 5 ~ 1, where 6 ~ IA2I. For pio and c^j given in Sect. 5 and Ho = 1, 
we obtain roughly Biq ~ 10~^ G. 

We have derived the ponderomotive force for circularly-polarized waves traveling along 
the magnetic field. For example, in papers by Denton et al. (2004, 2006) and Takahashi 
et al. (2006), toroidal Alfven waves were discussed. Therefore, it is worth to derive 
ponderomotive forces for other types of perturbations. The case of the strong nonlinearity 
when density perturbations are the same as background values deserves also to be examined. 

We have assumed that condition of quasineutrality is satisfied in the background and 
perturbation states. It is followed, for example, from the paper by Denton et al. (2006), 
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that equatorial peaking is observed for ions and absent for electrons. We think that a large 
electric field should arise in this case. 

It is obvious that for a detailed comparison of theoretical results with observational 
data, theoretical models should be adequate to real experimental conditions and observations 
and vice versa. It is important to be sure that a stationary symmetric theoretical picture 
considered here is relevant to real situations. The equatorial peaking in this picture is 
possible, if the wave action is symmetric and simultaneous from both ionospheric boundaries 
to the equator during some time for establishment of the stationary state and wave 
amplitudes should decrease from the ionosphere to the equator. In other cases, dynamical 
(nonstationary) processes can play a role. 
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7. Figures 
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Figure 1 




Fig. 1. — The dependence of the Earth's magnetic field B (in G) on the geomagnetic latitude 
(/9 in the meridional plane of the dayside magnetosphere for different L in the model by 
Antonova and Shabansky (1968) 
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Figure 2b 
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Fig. 2. — a) The dependence of the relative density perturbation 8 along the field line at 
L = 6 and Bio = 10^^ G. The curve 1 corresponds to Hq = 0, curve 2 - Hq = 0.5, curve 3 - 
Hq = 1.0. b) The same as in Fig. 2a at L = 10. 
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Figure 3b 
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Fig. 3. — a) The dependence of (5-distribution on x in the dipole {k = 0) and two-dipole 
{k = 13) geomagnetic field at L = 8 for Biq = 10^^ G and Hq = 0. b) The same as in Fig. 
3a for Ho = 1.0. 



-23- 



Figure 4a 
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Figure 4b 
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Fig. 4. — a) The dependence of the relative density perturbation 6 on x for different values of 
the wave amplitude Biq = (1, 3, 5) x 10^^ G at L = 10 and Hq = 0. The curve 1 corresponds 
to Bio = 10-^ G, curve 2 - 5io = 3 x 10"^ G, curve 3 - 5io = 5 x 10"^ G. b) The same as 
in Fig. 4a for Hq = 1.0. 
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Figure 5 
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Fig. 5. — The dependence of 5 on x at L = 10 and Biq = for Hq = 0.5. The curve 1 
corresponds to i7 = 0. The curve 2 - g^ff ^ 0. 



